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The  interaction  of  geomagnetic  ally  trapped  electrons  with  the 
neutral  atmosphere  entails  both  a diffusion  in  local  pitch  angle  (a) 
and  a loss  in  electron  energy  (E).  The  corresponding  transport 
coefficients  DOQ  and  (dE/dt)^  can  be  constructed  by  summing  the 
contributions  DQ^  and  (dE/dt)j  from  the  various  atmospheric 
constituents  (j).  The  subscript  v on  dE/dt  serves  only  to  emphasize 
that  the  corresponding  process  is  essentially  frictional,  as  distin- 
guished from  stochastic.  The  partial  transport  coefficients  D^^ 
and  (dE/dt)^  are  proportional  to  the  local  number  density  N.  of 
atmospheric  constituent  j but  are  independent  [e,  g.,  Walt,  1966] 
of  the  local  pitch  angle  (a).  Both  decrease  in  magnitude  as  the  energy 
(E)  of  the  incident  electron  increases. 

Radiation-belt  evolution  is  conventionally  described  [e.  g. , 
Haerendel,  1968]  by  means  of  a bounce- averaged  Fokker-Planck 
equation.  Thus,  it  proves  convenient  to  introduce  the  bounce- 
averaged  partial  transport  coefficients  D ^ ^ = <(Bq/B)  cos^s) 

and  <(dE/dt)j>,  where 

y = (Bq/B)1/Z  sin  o (1) 

is  the  sine  of  the  equatorial  pitch  angle  (oq)  and  Bq  is  the  minimum 
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value  of  B (the  magnetic-field  intensity)  along  the  field  line  (L).  The 
bounce  average  of  any  quantity  Q is  given  by 


= (n2/2w)^*(Q/v)  aec  or  da. 


where  v ia  the  speed  of  the  electron,  a is  the  coordinate  that 
measures  arc  length  along  J3,  and  is  the  full  bounce  period 

defined  by  setting  Q = <Q>  = 1 in  (2).  The  purpose  of  the  present 
note  is  to  derive  an  analytical  relationship  between  the  contributions 
of  a given  atmospheric  constituent  (j)  to  D and  <(dE/dt)^>,  i.  e. , to 
derive  an  analytical  relationship  between  <(Bg/B)N^  cos  «r>  and  <N^. 


DERIVATION 


The  essential  step  in  the  derivation  is  to  observe  that  (1 ) implies 
cos  a » [1  - y^(B/Bn)]*^.  Thus,  it  follows  from  (2)  that 


(2w/O2)<(B0/B)N.  cos  o>  = 


(Nj/v)(BQ/B)  ds 
[1-  y^B/Bg)]"1 


while 


y(2ir/fl2)  <N.>  = f - 


y (N./v)  ds 

[1  - y^B/Bg)]1 


Moreover,  the  right  member  of  (4)  is  equal  to  minus  the  derivative 
of  the  right  member  of  (3)  with  respect  to  y.  Thus,  it  follows  that 


2 /»! 

<(B0/B)N.  cos  «>  = n2(y)  I [u/n2(u)J  <N  > du, 

•'if 


since  cos  a s 0 for  y = 1.  This  is  the  desired  result.  The  practical 

Z 

■ignificance  of  (5)  is  that,  given  either  or  <(Bq/B)Nj  cos  a) 

aa  a function  of  y,  one  can  immediately  obtain  the  other  by  integra- 
tion or  differentiation  (respectively). 

It  can  be  shown  by  Jacobian  methods  [Haerendel,  1968]  that  the 
bounce-averaged  phase-space  density  f satisfies  the  Fokker-Planck 
equation  [cf.  Walt,  1966] 


n,  a I y°  af 


y 9y  ^ 9y 


77  « v v\w./l  (6) 


in  the  absence  of  range  straggling  (Dj.  = 0)  and  radial  diffusion 
(Dll  = 0),  where  Y is  the  ratio  of  relativistic  r -as  (m)  to  rest 
mass  (RIq).  The  present  work  has  shown  that  the  contributions  of  any 
atmospheric  constituent  (j)  to  and  <(dE/dt)^>  in  (6)  are  related 
by  virtue  of  (5). 


APPLICATIONS 

The  bounce  period  (2n/(22)  in  a dipolar  magnetic  field  (15) 
is  proportional,  at  fixed  energy  (E)  and  shell  parameter  (L),  to 
a function  called  T(y)  which  is  well  approximated  f Davidson,  1976] 
by  the  formula 


T(y)  * T(0)  - [T(0)  - T(l)]y' 


* 1.380173  - 0.639693  yJ/\ 
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Earlier  approximation!  of  T(y)  had  been  given  by  Hamlin  et  al. 

fl96l]  and  Lenchek  et  al.  [1961].  Although  T(y)  can  be  generated 

to  arbitrary  accuracy  by  a numerical  integration  of  (2)  for  Q = 

= 1,  the  use  of  (7)  in  (5)  should  make  the  required  integral  for 
2 

«Bq/B)N^  cos  a y much  easier  to  perform  in  most  instances. 

A particularly  simple  case  utilizing  (7)  is  that  of  a spatially 
uniform  atmospheric  density  N..  This  occurs  in  most  laboratory 
devices,  for  example.  (A  laboratory  device  that  simulates  the 
geometry  of  the  magnetosphere  is  commonly  called  a "terrella.  ") 

In  situations  such  as  this,  for  which  = Nj  (a  constant),  it  follow! 

from  (5 ) that 


<(Bq/B)N.  cos#>  = [Z(y)/T(y)]N., 


where 


Z(y)  = f u T(u)  du 

J \r 


a <1/2)(1  - y ) T(0) 

- (4/11)0  - y11/4)[T(0)  - T(l)] 


if  T(y)  is  approximated  by  (7).  The  present  approximation  yields 
Z(0)  a 0.  457471,  which  agrees  well  with  the  exact  result  [Schulz,  1974] 
that  Z(0)  = 16/35  (»0. 457143).  An  earlier  approximation  of  Z(y), 
which  yielded  Z(0)  a 0.455533,  had  been  derived  by  Schulz  [ 1974]  in 
a different  context.  It  follows  from  (6)  that  the  Fokker -Planck 
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equation  for  scattering  by  a uniform  atmosphere  is  of  the  form 


D 3 df 

— ——  y Z(y)  — 
y T(y)  dy  ay 


V v 9E 


f . (10) 


where  and  <(dE/dt)^>  are  independent  of  y and  of  L. 

In  the  real  magnetosphere,  of  course,  one  does  not  encounter 

anything  like  a uniform  atmosphere.  However,  various  investigators 

have  numerically  implemented  the  bounce  averages  indicated  in  (3) 

and  (4).  The  dashed  and  solid  curves  [Walt,  1966]  shown  in  Figure 

2 3 3 2 

1 represent  (respectively)  the  quantities  (2  Y v /c  x ) D and 
(-v/mQc3)  <(dE/dt)w>,  where  x is  the  cosine  of  the  equatorial  pitch 
angle  (hence  x + y = 1)  and  c is  the  speed  of  light.  These  par- 
ticular forms  serve  to  factor  out  the  energy  dependence  of  the  transport 
coefficients,  except  for  some  slowly  varying  logarithms  which  are 
evaluated  at  E ~ 1 MeV. 

Except  for  the  fact  that  the  neutral  atmosphere  has  several 
constituents  (j),  one  should  expect  to  find 


(v/m0O<(dE/dt)v> 


yT(y> 


-I 

dy  ( 


(1  - yc)  T(y) 


' 2Y2v3  1 ) 

.~77d”\  ) ' 


where  A is  a dimensionless  constant  that  depends  in  detail  upon  the 
ratio  between  the  inelastic  and  elastic  cross  sections  of  atmospheric 
molecules  or  atoms  for  incident  electrons.  The  "data"  points  (filled 


circles)  in  Figure  1 serve  as  an  illustration  of  this  idea.  The  square- 
bracketed  factor  in  (11)  has  been  differentiated  numerically  at  L = 1.  185 
and  treated  as  a constant  (for  x < 0.86)  at  L = 1.9  for  this  purpose. 

The  resulting  "data"  points  (normalized  by  evaluating  A to  force 
agreement  at  x = 0)  follow  the  solid  curves  rather  well.  One  would 
not  have  expected  perfect  agreement,  since  the  densities  (bT)  of  the 
various  atmospheric  constituents  (j)  have  different  altitude  profiles 
and  contribute  to  D and  <(dE/dt)..>  with  different  weights  (essenti- 

yy  w 

ally  as  Z2  to  and  as  Zj  to  ((dE/dt)^>,  where  Zj  is  the  nuclear 

charge).  Thus,  the  transition  from  hydrogen  dominance  to  oxygen 
dominance  along  field  lines  on  the  L = 1.9  drift  shell  may  help  to 
account  for  the  departure  of  the  "data"  points  from  the  solid  curve 
at  large  x. 

Bounce-averaged  atmospheric  densities  have  been  compiled 

by  several  investigators  [e.  g.,  Cornwall  et  al. , 1965]  in  the  course  of 
studies  on  radiation-belt  protons  (for  which  atmospheric  pitch-angle 
diffusion  is  an  unimportant  process),.  The  present  result,  as  summari- 
zed by  (5),  enables  one  to  extract  the  weighted  bounce  average  relevant 
to  electron  pitch-angle  diffusion  directly  from  such  compilations  of 

For  example,  if  one  has  found  that  £nNjj^y  n for  some 

range  of  L values,  then  it  follows  from  (5)  and  (7)  that 


<(Bq/B)  N.  cos2o> 


£ 


N' 


f(j) 

n 

T(y) 


1 - 


•n 


2.75 -n 


(2  - n) 


0. 53939  (11  - 4n) 


(12) 
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Terms  corresponding  to  n = 2 and  n=ll/4,  if  present  in  (12), 
reduce  to  logarithmic  form  upon  expansion.  Applications  such  as 
this  may  enable  the  present  discovery  (which  seems  just  mildly 
curious  at  first  sight)  to  be  useful  in  practice  as  well  as  in  theory. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 


Acrophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer. reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 


Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos* 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 


Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propsgstion  phenomena,  including  plasma  electromsgnetics;  quantum  electronics 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 


Materials  sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor 
rosion  and  fatigue-induced  fractures  in  structural  metals. 


Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion irom  tKe  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms , and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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